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ABSTRACT 


The contribution of electron exchange to intermolecular 
potentials is studied using the operator form of perturbation theory, 
when the non-interacting system of atoms or molecules is degenerate . 

A general, symmetry adapted degenerate perturbation theory is 
developed, including the almost-degeneracy case* Perturbation 
expansions are given up to second order in the energy* Differential 
equations for the perturbed wavefunctions are obtained, as well as 
their associated variational principles. 
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I. Introduction 

When dealing with the interaction of atoms or molecules, fcr 
example in collision processes or inside crystals, the wavefunc tions 
and energies for the non interacting system are sometimes accurately 
known, in which case it is convenient to use perturbation theory to 
obtain the corrections due to the interaction. Furthermore, the 
new wavefunctions are required to have a symmetry corresponding to 
an irreducible representations of the group of the full hamiltonian. 
This symmetry is usually imposed by applying a group theoretical 
projection operator , for the appropriate representation, to the 
set of unperturbed states. But a set of linearly independent 
unperturbed wavefunctions may become linearly dependent after being 
projected, so that the development of a perturbation theory in 
terms of the eigenstates for the non interacting system requires 
special attention. 

This problem is present in the perturbation theory of inter- 
molecular potentials in the intermediate range between long distances, 
where a multipole expansion is used and exchange effects neglected, 
and short distances, where chemical binding is found and the exchange 
effects predominate. Here the operator JL antisymmetrizes the 
total wave function with respect to electron exchange. 

The perturbation theory of intermolecular potentials including 

exchange was initiated by Eisenschitz and London'*' and continued 

2 -8 

since then by several authors . The use of linearly dependent basis 
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sets has led sometimes to inconsistencies in the literature, which 

9-12 

have been avoided more recently by using linearly independent, 

usually or thonormalized, sets. 

13 14 

Another fruitful approach y has been based on operator 

15 16 

perturbation theory, 9 which helps to isolate the difficulties derived 
from the symmetry problem^’. In this contribution we want to use the 
operator approach to study the problem when the non interacting system 
is degenerate 0 For example, degeneracy will exist if the non 
interacting atoms or molecules have orbital degeneracy or their 
spins coupled to a multiplicity smaller than the highest:., We shall 
keep these examples in mind but develop the formalism in general* 

In section II we give the basic equations for symmetry adapted, 
degenerate perturbation theory and consider also the almost- 
degeneracy case, Section III includes the perturbation expansions 
when the degeneracy is completely removed in first order and when 
it persists to first order and is removed to second order. Section IV 
gives the corresponding differential equations and discusses their 
associated variational principles. 

II. Symmetry adapted, degenerate perturbation theory 

Let us consider a system of atoms or molecules whose energy 
in the absence of interactions is g tim.es degenerate, with the 

of linearly independent, non-symme tr ized wave functions will form, 
together with the set of excited states j k > 8 , which 


§ 




This set 


unperturbed states given by ^ C|)^ ^ ^ 
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verifies <C k l ' ^=*0 * a complete basis to describe the 

non-interacting system. 

We want to obtain a perturbative solution for the interacting 
system in terms of this complete basis^ including the symmetry 
properties of the whole system* Indicating with the group 

theoretical projection operator which selects the wavefunction 
components belonging to the irreducible representation of interest, 
it will in general be a product of a point group projection operator 
times the antisymmetry projection operator. The set of symmetry 
adapted states n0t necessar ^^y linearly independent 

and could lead to ill behaved expansions* In order to avoid this 
difficulty we introduce the g x g overlap matrix (or Gram’s 
matrix) A for the set^ 


a = C I q> CJ >3 , (1 

which is hermitian, and diagonalize it by means of the unitary 
transformation S to obtain the diagonal matrix y 


S+A 


S = % 

ro — /-v/ 


( 2 ) 


The eigenvalues (V- of A are definite positive, i.e. ^ O 

\ S\J “i / 

and we indicate the corresponding eigenvectors^ the i-th columns 


of S y by S. , so that 

a; 1 


A Sy = <S j 

f>jJL V 


( 3 ) 



Selecting the h ,$ g eigenvectors with eigenvalues different from 
zero we form the states 



where h^p^l and Cp o is a row matrix with elements 
These states form a set of symmetry adapted, linearly independent 
wavefunctions describing the non-interacting, system of atoms or 
molecules* In effect, clearly , and 


S* 5 31 <cf ol | ofy> 
= ( <v Wj ) 2 C S+A S )^_ 


so that they are orthonormalized . Notice that this result is 

obtained without imposing the orthonormality of the set * 

We shall use the previous results to set up the eigenvalue 

16 

problem using the wave and reaction operator formalism , indicating 
with capital greek letters the symmetrized states, and introduciijg 
the projection operators 




( 6 ) 



and 


9 = . ( 7 ) 

In Schrodinger equation ch-e)tjt» 0 7 = , 

we have and Cf-M! = O so that only the 

projection with ^ will give a right hand side different from zero 
and we can write., in the usual way, 

&CH-E + C£CH-'E)'P<Ur = 0 J (8. a) 


c ?C«~E)fc>Tf +T O . 


From (8„b) we obtain 


Tw TH <S>^£~ , T =9 [W d-jl )y>to +9 CE -H )?]- | p 

(<0 

where ^ and insure that T is well defined and do not 

affect the results, and replacing (9) in (8 „a), 

to W-E)Wto-tjJ-=o 7 W = Co+TH . ( 10 ) 


Introducing the notation 


^ <^ p |^r> = Z- 


(ii) 
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and projecting the first Eq. (10) on 
secular equation 






we obtain the 


IHwlx^p-^ , (12) 


discussed in detail by Lowdin for ^ = 1 . Here the index ^ ^ 1 
numbers the energies that split from due to the interaction, 

and we have used the result K ( W I ^ 

Once the coefficients are known the states corresponding 

to are obtained from (9), 

a 


fo^+'7^? 1 = 2- a+Tfi )^ 0( , 


(13) 


Equations (12) and (13) may be rewritten explicitly in terms of the 
unsymmetrized states 


using the decomposition H = H Q + V , to obtain the secular 
determinant 


C< s ) 



7 


For i = p = 1 we put ^ oA - Cfo ? ‘f© Ijfl I j 

and re-obtain Van der Avoird*s result for the non-degenerate case. 
For h = g, i.e., when the set Is linearly independent, 

the full JS unitary transformation appears in (15), which can then 
be eliminated to obtain, in matrix notation. 


<*«* [<^|v^|q> 6 >+ ^ivn/l^-CE-e^A 1 

rO r* c*J ^ ~ 



If the matrix elements of second and higher order in V are neglected 

in (16) we are left with the valence bond secular equation usually 

obtained when the interacting system wavefunction is written as a 

projected linear combination of Slater determinants constructed from 

18 

the products of isolated atom wavefunctions 

It is seen from eq. (12) that we can always choose 


V" 

C . c _ Kl 5 
t If 1f“ l 5 !'! 


(17) 


and using (13) 


-t. |^> or) 

2quation (12) together with will then determine 
the coefficients } £, ^4 3 uniquely or alternatively, 
thoosing 


N will determine the norm of It is also 

1 * 
from Eq. (18) that for fixed c n values the states 

4 - 


seen 


are 



orthogonal within an error of second order in V . In the following 
section we shall find it convenient to use the condition 


21 c^ + c a = i 
t *t* IP 


(19) 


,<o) 


where the are obtained from the lowest order solution of 

Eq. (12). 

In some cases the unperturbed state of interest will not be 
degenerate, but its energy will be close to other unperturbed energy 
levels, their difference being small compared with the magnitude of 
the interaction. This situation may arise when dealing, for example, 
with spin-orbit effects. These almost degenerate cases can be 
treated in the way described above^"? Let us indicate with Ef ^ and 
^ ^ i <L , the closely spaced energies and their wave- 

functions and let ^ , k>g , indicate the other unperturbed 
states. We introduce the operator ^ defined by 


A- ZAiKfoXtM 


(20a) 


3 j-' j A » 


Cc>) 

C 


(20b) 
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We find then that 


•Mo; 


: C*) 

'A\t T 0 ^ 


% 




? 


( 22 ) 


so that is 8 times degenerate and is a g- order 

degenerate set. The previous and following results may then be 
applied to these problems by simply replacing Eq , and V 


by E 


Co) 

Air 


, Hq and V respectively. 


III. Perturbation expansion in the operator formalism 

To obtain the eigenfunction and eigenvalue E of the 

^ % 

interacting system of atoms or molecules to increasing order in the 

potential V , we require a perturbation expansion for T . This 

14 

we shall do, following Van der Avoird's treatment , by noticing that 


= ~-A- Vjft 


(i) 


where -A. is a normal, idempotent operator defined by 


TV = -i - ^ , 


( 2 ) 


which verifies 


, and where 


V’ = V— ) . 


( 3 ) 
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From (1) and the defining equation for T , 

T OE-+OT= (^ } -VIo)T-A\/'T=?. (4) 

(5) 

(6) 

(7) 

Using the relation 

(si 

and iterating in Eq. (7) it is found that, provided the iteration 
expansion converges, T = T ' . The hermitian character of T may be 
explicitly shown by replacing T in the right hand side of (7) by 
T + , so that 

4-Tr^v't v!a + c 0 T 

=?i? 0 T+ '? -A. ^ C V'i)+ ^ v0At'P o ?+ 'tv‘A~ <?JP 

CO 


Introducing the unsymmetrized resolvent operator 

l?o- f- \tp *cf k | 

satisfying 

and operating with f 3 >t Ro to the left of (4)^, we get 

+?r„av' t 
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In this form T may be expanded as a sum of hermitian terms , which 
may be obtained in the basis set of unsymmetrized states for the 
non-interacting system. 

To calculate and we must solve by iteration the 

\ 1 

equations 

\ +<^ r !VTV Vp%° 

and 


= TVlc 1 = TV/4> . < 10b > 

1 If 7%. 


The procedure to follow will depend on the degeneracy being completely 
or partially removed, or persistent in first order, or second order, 
etc. We shall only consider two cases: 


I. The degeneracy is completely removed in first order. 

II. It is persistent to first order but completely removed to second 
order . 

Other cases may be dealt with in similar ways. 

We start with case I. To first order we get from Eq. (10a) 

t (< K' v f ) ^ = ° (u) 

which gives h different eigenenergies E ^ and, choosing 


12 




( 12 ) 


the unsymraetrized states 

Vj’.iJcpw 

(13) 

and using 

TVc^ = j\ v (V-eJ ) ^ ( u) 




ZL 

% “ f> 

Replacing 


,o) 


fo 


'If l°f 




and 


in (10b) we get 


the result is 

(V-E f ' ; )CP^ • < 15 > 

-T 

To obtain the coefficients & a required in fo^° , as 

' ^ (2) ** 

well as the second order energy E A j we proceed to the second 

+ 

order approximation for (10a ) > given by 


<+* IV3>^’>+Z (<^ r \s/fi 


_ e 


(2) 

1 
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y- 

Multiplying by ° , summing over p' and using Eq. (11) and 


(13), 




f?> + <*>**£ 


— U) 

~ fc <l 5 ^ = '° 


C-ITJ 


so that for 0 I =r 


X=% > 


«?JV- 


(18) 


To find we write 

'4 “45 1 

and replacing the second expression in Eq. (17) we get, for ) 


or c 0, ,«Z-v-« (19) 

n‘ y\%< i'f 


^ = ^f E 4' ) & ji c v-e«) j<f 


f<0 


'1 


> 

(2oJ 


fl) 


-I* 


r v \ ) />\V 

^ =0 , resulting from 22 C J 

* ^ 

obtain from Eq. (10b), 


, define . We also 
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rj? v cf 

+ ?R„Ai [CV-^’)jA4 J! CV-E^)]^ R B ^Cv-w“ , j<pJ ) 

C*0 


This procedure may be continued to obtain higher order results. For 
example, the third order energy will be given by 




i CCV-E ( '')i +fi A+R = jH (V-^)lCfp 


(JL%) 


In case II we would get, from the solution of equation (11), 
= , for h^^'^d , and the 

would be yet undetermined up to an unitary transformation. We go 


s tates 


then to the second order approximation to Eq. (10a) and write 


+- 


f ' V - - ^° s ff 3 - o 


02*) 
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Putting 


c <0 ' 0) = Z d io) c lo) 

^ X' W ? 


( 24 ) 


multiplying Eq. (23) by /C l . ( 


after using Eq. (11), 


(c>)“^ 

4 1 p‘ 


; and summing over p 1 we obtain; 


X ( I V ?fV? v I f ^ l^T° (25) 


4 tv n* 


fO,£) 

which gives us h different second order energies E , 

4 


and 


f ' = Z Jr = Z . cp^> 
1 4 f it> ’ 6 |“ 41 41 14 ' 

i <f 


(26) 


( 0 ) 0 ) 

4 


We also obtain from Eq. (10b); 


^ 1 T 


(27) 


To find /C C °' 0 and E £ °»^ 
in case I; and so on. 


we would go to the next order; as done 


This completes our discussion of the perturbative solutions » 

To compute second order energies we would need to perform sums over 
the excited states of Hq ; included in Rq . When dealing with 

corrections to ground state energies; it will be useful in some 

, . ’ . 16 
cases to use the inequality 


16 


o > i*« > ( )£ if b y% i (28) 

3 fe-g-M 


from which the following bounds are obtained for 


nW 


) 


0>^1 


( 2 . ) 

% 


> | cv/-e" 4 ')i 

* ^CV-E«»)IC|>W> . Ci1 ) 


The Unsold approximation may be considered as an interpolation 
between these two bounds. 


IV. Differential Equations and Variational Principles. 

In connection with practical applications ^ it is frequently 

useful to work with the differential equations for each perturbation 

20 

order and their associated variational principles 

To obtain a differential equation we shall look for a wavefunction 
belonging to the total (unsymmetrized) space^ such that 




to 


From the relation 


v't 


(2. a) 


T w =7 D* a-j )a[MS) +7 (*>'- u>)T Y { P 


(2.b) 
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and (III. 10b) we get, introducing $ -A— 7 


- T Co Vcf^+ + Cg^ , < 3 > 

and operating to the left with < " _ P -H 0 ) ) 

?(E«-« 5 )^=?v^ + ?V'}A+^ . <4, 

Rather than working with this equation we may try to solve the 
"related" equation 


CE?- Wo ) ^ =Tv^ + 7 v'a+ 

and, since this would also solve Eq. (4), we can write 

. The normalization condition for requires 

that < be finite. Otherwise we are free to 

choose, to solve Eq. (5), the boundary condition Cjp^ J ^ o^)> 

9 > 4 . It is then seen that the formal solution 

3 ^ " 


(5) 


^ = T?.7i/fy + At ^ 

= <w- ) _l i? t dP v q>. 


gives, after expanding and rearranging terms, a solution for Tya 
identical to the one obtained in Section III. Equation (5) together 
with the adjoint related equation 
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( £?- Yfa J?- V ^ + yv vijT^ <■ 

may be obtained from the variation function 

(< ^? lTV 1 V 


+ I f % >' < \ I ^>) 

* Cf) 

where < Cfoi I ^ > = <* %X. I % £ >3^ > ' L ^ * > 

and the variational principle 


s V s %3= s J 4 /s ^^ = o 


(9) 


It is seen from Eq. (8) that in order to get an< ^ ^ 

necessary to solve first the secular equation (Ill.lOa)^ to obtain 


E 4 and 




1 


This may be done using again a perturbation 


-^(0 j n>^ 
treatment. Once fc ^ and L r 


are known from the solution of 

% """ • * , () 

Eq. (III. 11) one would obtain from the first order 

' 1 

approximation to Eq. (5), 




(10) 


for case I of the previous section. A variational expression for 

o=> 

may also be found expanding JT- — 21 ^ and 

> + 




*1 


= o * 
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. . , TW T CC' 

noticing that J n — J , 


X ^ ' — O • Since Eq. (10) does not 
include the non-hermitian operator 

-r-(2) 

and 0 ,, may be written 

X 


it will be 


y 4 C fip = < Xj° 1 ^4 > + iv/? 


— <v (l) 




OH] 


where we have used ^ ^ ^ = <1 

91 * 

Variational 

principles for higher orders and for case II may be discussed along 


^ 1 1 ^ 

symmetry adapted Hylleraas variational principle 


Equation (11) is a 
21 * 


Putting 

found that t £«> 3 > J !' C j 

so that (pleads to a minimum variational problem for JE 


in this equation it is 

(I) —t -f — 

J= £ A ? 

fe) 


^ T 


4 


1 


the same lines 

Variational principles of the Hylleraas type, like Eq. (11), 
are useful provided the solution to the zeroth order problem is 
accurately known. When this is not so it will be usually more 
convenient to use the full variational principle. It is also seen 
that charge transfer effects may be studied by incorporating in Kjj, CM) 
the corresponding charge transfer states. Nevertheless^ that may 
complicate the computation of the last term in Eq. (11) ^ making 
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again more convenient to use the full variational principle. The 
best computational method to follow should be chosen in accordance 
with each particular problem. 

V, Discussion 

In the previous treatment we have assumed that the solution for 
the non-interacting system of atoms or molecules is accurately known. 
This will certainly be the case when dealing with interacting 
hydrogen atoms. The equations obtained for the energies and wave- 
functions may be used to describe a pair of hydrogen atoms in a 
symmetry state corresponding to excited, i.e. orbitally degenerate, 
non-interacting hydrogens, or to describe more than two hydrogen 
atoms in a state with multiplicity smaller than the highest. For 
other atoms the usefulness of the perturbation approach with depend 
on the accuracy of the unperturbed states. 

As it is usually the case with perturbation theoretical 
treatments, we have assumed that the iteration solution of the 
equations converges, in which case the energy has been shown to be 
real to all orders. With respect to the meaning of "order" in 
symmetrized perturbation theory it may be useful to take into 
account that the zero order problem is recovered by taking V'— > O 
and ^ 4 . For example, when dealing with Eq. (IV. 5) we can write 


where in the right hand side we apparently have a zeroth order term. 
Never the less, for the operator commutes with 

and the term disappears in the unperturbed problem. Consequently, 
the whole expression must be considered as a first order one. 

From the knowledge of accurate to a certain order it is 

+ 

possible to compute expectation values of properties of the inter- 
acting system to the same order. The applicability of Dalgarno ! s 

interchange theorem for expectation values is restricted as in 

20 

unsymmetrized degenerate perturbation theory . This restriction 
might have been expected, since our results reduce to those of 
unsymmetrized degenerate perturbation theory f or j4 = 1 
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